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CN . We point out the existence of new effects of global spacetime expansion on 

local binary systems. In addition to a possible change of orbital size, there is a 

contribution to the precession of elliptic orbits, to be added to the well-known 

general relativistic effect in static spacetimes, and the eccentricity can change. 

Our model calculations are done using geodesies in a McVittie metric, repre- 
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CN ■ senting a localized system in an asymptotically Robertson- Walker spacetime; 

J^ , we give a few numerical estimates for that case, and briefly comment on ways 

in which the model should be improved. 
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The issue of whether the global cosinological expansion affects local gravitating systems, 
such as planetary systems and galaxies, has been studied for a long time (for recent work, 
and discussions of previous results, see e.g. Refs. |]l]-§]). The specific questions asked in 
the literature focus mostly on the extent to which local systems expand, and on the form 
and magnitude of the corrections to the effective forces felt by orbiting test bodies; while 
the approaches used vary, sometimes in conceptually important ways, the consensus is that 
there is an effect in principle, but in practice it is exceedingly, undetectably small. To a first 
approximation, a reasonable, physically motivated point of view, as expressed by Misner, 
Thorne, and Wheeler [Q, for the case of galaxies, is that they are like rigid pennies attached 
to the expanding balloon representing the universe, and they do not themselves expand. 

While we agree with this general conclusion, at least in terms of currently feasible ob- 
servations, an increase in orbit size is not the only effect the cosmological expansion can 
have on local systems, nor the only cumulative one. In the search for such effects described 
in this paper, we use as a model spacetime an actual solution of Einstein's equation which 
represents a gravitating body in an asymptotically expanding universe; restrict our atten- 
tion to nearly Newtonian orbits in the weak-gravity, slow-expansion approximation; and 
consider the time variation of the parameters characterizing the corresponding Keplerian 
ellipses, their size, perihelion angle, and eccentricity. As we will show, this approach leads 
us to predict changes in all of those parameters as a consequence of global expansion. 

Before we proceed, however, we will make a few comments on the systems involved. It 
is well known that Keplerian elliptical orbits in Newtonian gravity do not precess (the 1/r 
potential is one of two that yield closed orbits, together with the harmonic oscillator, as 
Bertrand's theorem states) but the corresponding ones in general relativity do, while the 
size and eccentricity of the general relativistic orbits are constant. These results hold for the 
case of a single, spherically symmetric, isolated center of attraction, with a test body and no 
other matter around it. Corrections due to other orbiting bodies can be calculated, and they 
may affect all parameters of the orbits, but as far as the global expansion is concerned, under 
the above conditions there can be no local effects due to cosmological expansion in general 
relativity, regardless of the rate at which the rest of the universe expands, because Birkhoff 's 
theorem then concludes that locally the spacetime must be the static, Schwarzschild solution 
(Einstein and Straus also showed this explicitly in one class of models, without using 
Birkhoff 's theorem). The cosmological expansion can be felt by a local system only if (i) 
The situation is not spherically symmetric, and/or (ii) There is matter (some gas, dust, or 
dark matter, say) within the local system, or a cosmological constant. This is one of the 
many ways in which one can see that general relativity obeys the spirit of Mach's principle, 
but only up to a point (see, e.g., Ref. [0]); the global behavior of distant objects can affect 
the orbital dynamics, but the dynamical equations are local. 

In a realistic model of a local system, both deviations from spherical symmetry and extra 
matter will have an effect; in the solar system, for example, the presence of other stars, solar 
multipoles such as oblateness, minor objects, interplanetary matter and/or the solar wind 
and radiation, are ways in which this happens. Studying the way the global expansion 
affects the local system then means: (i) Understanding how the dynamics of surrounding 
non-symmetric matter and/or the local internal environment of the system depends on the 
global behavior of the universe, and (ii) Calculating the effect of the local environment on the 
dynamics of the system of interest. Notice that, for a local system that is considered as part 



of a larger one, such as a planetary system inside a galaxy, the first part amounts to assuming 
that a similar problem has been solved one level up; for example, the spacetime metric for a 
planetary system should not be asymptotic to a Robert son- Walker metric expanding at the 
Hubble rate, but to a metric that is appropriate inside the galaxy, expanding, if at all, at a 
rate previously found for that system. The problem thus becomes somewhat involved, but 
it seems to us that this is the only way to avoid sneaking in an assumption about the local 
environment's expansion as part of the setup, which has occasionally been done. On the 
other hand, once it has been formulated in this way, the problem is broken down into parts 
that can be separately investigated, and labelled by the structure scale under consideration 
(galactic, planetary, ...), mechanism by which the effect may be transmitted, and type of 
effect. 

Here, as a start, we will limit ourselves to considering a simplified model to show the 
effects that may arise. Spherical symmetry will be preserved, and the mechanism responsible 
for transmitting environment effects will be the presence of matter in the local system. In 
our numerical estimates we will use orbital parameters relevant for a couple of known planets 
and galaxies, and the Hubble constant to obtain the environment expansion rate; much more 
work will be needed on several aspects of this problem, including the type of matter and 
spacetime model used and the role of anisotropy, before we can claim to have a reasonable 
understanding of it along the lines described above. 

While spherically symmetric, vacuum, asymptotically flat spacetimes and homogeneous, 
isotropic cosmological ones with fluid matter or a cosmological constant can be easily 
treated in general relativity and give rise, respectively, to the Schwarzschild solution and 
the Friedmann- Robert son- Walker or de Sitter spacetimes, solutions representing an isolated 
massive object embedded in an expanding universe are more difficult to get, and fully ex- 
phcit forms are not usuafly given (see, for example, the general discussion in Ref. [^). The 
known metrics representing isolated bodies in asymptotically FRW spacetimes are mainly 
of three types: McVittie M, Einstein-Straus H, and Ferraris- Francaviglia-Spallicci [|TU|. In 



this paper, we will use McVittie metrics, although an explicit solution is given in Ref. |[10| 
where a Schwarzschild and a FRW solution are matched smoothly across an extended region 
of spacetime represented by a transition metric, where the orbit of a test particle may also 
be analysed. 

The (asymptotically) spatially fiat McVittie model with k = 0,m the isotropic coordinate 
form used by Hogan [|ll|], has line element 

ds^ = - ^7;/^ dt^ + e^(*) 1 + dr^ + rMfi^ . 1 

Here, the constant M is interpreted as the "mass at the singularity," and /5(t) represents the 
(asymptotic) expansion rate of the universe. We will not give an explicit expression for /5(t), 
since it is related to the specific equation of state of the fiuid, but will think of it instead in 
terms of an expansion /5(t) = l3o + Poit — to) + 0{(t — to)^), around the present time t^, say, 
where the first few coefficients can be fitted phenomenologically. 

The above line element has the advantage that, as r ^ cx), the t = const hypersurfaces 
become the surfaces of homogeneity of a Robert son- Walker model. For our purposes however, 
we find it more convenient to work with a radial coordinate with a geometrical meaning tied 
to the area of the corresponding 2-sphere, and use R := re'^/^(l + G'Me~^/^/2r)^ in the 



interval R G (2M, oo) 0. We can then rewrite the hne element as 
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cMt^ + ^^ -^ cdt dR + ^ + R^ dQ^ , (2) 



where /(-R) := 1 — 2GM/{c^R), and we have restored all c's and G's. 

The use of the McVittie metric as a model for studying local systems in FRW universes 
has been the subject of criticism [0,|1^], but in a series of papers |1^,^,|T5| Nolan made a 



good argument for the k = model by studying its global properties in detail. Nolan also 
showed that, in an appropriate sense, those solutions are the unique ones representing an 
isolated mass surrounded by a shear-free perfect fluid in a spatially fiat FRW universe; the 
shear-free condition leads to special properties that we would not expect a black hole to 
have 0] , but the metric does provide us with a viable explicit model for our purposes. 

Because of the spherical symmetry, we consider as usual our orbits to lie in the sin 6 = 1 
plane, and drop 6 from the whole treatment from now on. For a particle of mass m moving 
in the metric (§), one way of deriving the equations of motion is to use the Hamiltonian 
approach. The canonical momenta obtained from the action S[x] = {m/2) ^ dr g^i^x^i." are 
j9^ = mg^ijX^ , where the overdot on a dynamical variable, t, i?, or 0, denotes a derivative 
with respect to proper time r along a trajectory (whereas /J just indicates the derivative of 
the known function (3{t) with respect to its argument). 

The test particle Hamiltonian Ti = {2m)^^g^'^p^pjj in this case becomes 
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supplemented by the constraint that g^^p^Pi, = —{mcY, or 7i = — | mc^. 

Angular momentum is conserved, p^ = L, and as is common in celestial mechanics orbit 
problems, we reexpress R in terms oi u := 1/R in the equations of motion. If we perturba- 
tively expand those equations in powers of 1/c and f3, zeroth order terms correspond to a 
particle in Newtonian gravitation; neglecting all /? terms corresponds to doing calculations 
in a static spacetime, and we recover the well-known Schwarzschild results; the expansion 
effects we are interested in arise when leading order terms in f3 are kept. 

The radial equation can be recast into an orbit equation in terms of m(0), 

„ _ GMrn^ 3 GMu'^ /3'^m'^ CM /3m , 
"" '^''~~L^ + ^ 4Z^~ 2c^Lu "" ■ 

If we look for a solution of the form u(0) = uq + ^{(p), where uq = GMw? / L"^ corresponds 
to a Newtonian circular orbit and ^{(j)) is a small perturbation (see, e.g., Refs. [|T6|jr7| ), the 
linearized equation for the deviation ,^(0) is 

, _ "iGMul _ $^m^ _ ( QGMuq 3/5W\ GM$m 

It should be noted that for a general McVittie metric, /3 = $(t{(j))), so, to simplify the 
solution of the equation, we will now limit ourselves to considering the case /3 = 0; in terms 



of comoving coordinates as in Eq. (|T]), this is the hnear expansion case, although the metric 
for this model has a timelike Killing vector field, as can be seen from the fact that the line 
element (Q) becomes time-independent (in fact, it becomes equivalent to a Schwarzschild- 
de Sitter solution Q]). Equation (|) is then of the form ^" = fc — ^^ ^ + fd(,', which can be 
viewed as that of an oscillator with a constant force /c, modified by a /3 term which affects 
the equilibrium displacement of the oscillator and thus the geodesic orbit size, but it does 
not lead to orbit expansion (when /3 = 0); the frequency a; is also modified by a /3 term, 
which contributes to the orbital precession, and, interestingly, there is a damping term /d, 
which leads to a change in orbital eccentricity. The equation admits a general solution of 
the form ^(0) = C + Se'^'^ cos^Acj) — (po), depending on arbitrary constants £ and (po, and 
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from which, defining Uq = uq + C and Eq = S/{uo + C), we finally get for the inverse radial 
coordinate 

u{ct>) = u'o[l + eoe^^cos{A^-^o)]. (4) 

The angle a by which the orbit (^ precesses during each revolution can be found from 
the fact that the change A(j) such that the argument of the cosine increases by 27r is 27r + a = 
271/ A] i.e., for ^ ^ 1, o- ^ 27r (1 - A), or 

_ GttGM Svr $y{l - ey 

""^ a(l-£2)c2 +T GM ' ^^ 

where e = Eq e^"^ is the orbital eccentricity, and we have used the fact that for a Keplerian 
ellipse 1/uq = a (1— e^), with a the semi-major axis. The first term in (^) is the Schwarzschild 
contribution, o"o, as expected; the second term is the expansion-induced one, and has the 
same sign as the first term. The eccentricity also changes, by a fractional amount per 
revolution Ae/e = e^'^^ — 1, or 

Ae nGMpom tt/Jq 



2c'^Lun 2c2 



GMa(l-£2). (6) 



Equations (^ and (^ are our main results; they show that, even in a model with no orbit 
expansion (this may just be a feature of the /? = McVittie models used), global expansion 
has an effect on periastron precession and eccentricity. 

We can get a first crude estimate of the magnitude of these effects by evaluating a and 
Ae/e for a few known systems. For the solar planets Mercury and Pluto, even using a Hubble 
parameter value Hq around 75 km/s/Mpc for Pq, which is certainly a vast overestimate, 
one obtains for the relative size of the expansion-induced precession term in (Kl), values 



(jp/uQ Ri 3 X 10~^^ and 7 x 10~^^, respectively, and for the relative eccentricity change in 
(I), Ae/e ?» —1.5 x 10~^^/rev and —1.5 x 10~^°/rev, respectively. From the M and a 
dependence of these quantities we see that none of the strong gravitational field situations 



in which other general relativistic effects are studied, such as the binary pulsar ||T8[ or stellar 
motion around the galactic center |T^ can give usably large values, and we conclude that 
all expansion effects are negligible for the astrophysics of stellar systems. 

What about gravitationally bound pairs of galaxies? After all, as one might expect from 
cosmological effects (and in contrast to the Schwarzschild ones), both a^/aQ and Ae/e grow 
with orbit size a. We can get an idea if we apply equations (^ and (^ to the Large Magellanic 
Cloud, at a ^ 50 kpc from the Milky Way, with M k, 10^^ Mq, which gives ap/aQ ^ 0.34 
and Ae/e ~ —8 x 10^^^/rev, both significantly larger than the previous values, by many 
orders of magnitude in the case of ap/ao, in a regime where the use of Hq for /5, although 
probably still not appropriate, comes closer to being realistic. 

So, if some of the predicted numbers are not necessarily small, are our expansion effects 
cosmologically relevant? The problem here lies with the time scales involved in galactic 
dynamics, which make any direct observation of time variations impossible. Still, the fact 
that the Hubble distance-redshift relation shows local deviations on small cosmological scales 
is an indication that on those scales local interactions and global expansion effects coexist. 
Similar ideas have already motivated, e.g., work on the effect of global expansion on the 



formation and evolution of clusters of galaxies [20|, and more recently on the effect of 



inhomogeneities on the overall expansion [^ or of a cosmological constant on local dynamics 
2^ , and one of the potential benefits of this line of research is the possibility of gaining a 
new handle on dark matter. The present approach can be seen as a model for the onset of 
this situation with simple, binary systems, and as we have seen it can provide predictions for 
such systems. The model can be improved to include more realistic sources, diffuse matter 



and anisotropy, as well as possible multipole moments of the orbiting mass |2^ . It may then 
be possible to identify statistical effects which can be measured by observing large numbers 
of galaxy pairs, or early universe effects which left and imprint on later evolution. 
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